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We investigate how a magnetic field induces one-dimensional edge channels when the two- 
dimensional surface states of three-dimensional topological insulators become gapped. The Hall 
effect, measured by contacting those channels, remains quantized even in situations, where the 0- 
term in the bulk and the associated surface Hall conductivities, a^y, are not quantized due to the 
breaking of time-reversal symmetry. The quantization arises as the 0-term changes by ±2nn along 
a loop around n edge channels. Model calculations show how an interplay of orbital and Zeeman 
effects leads to quantum Hall transitions, where channels get redistributed along the edges of the 
crystal. The network of edges opens new possibilities to investigate the coupling of edge channels. 



In topological insulators (TIs) the topological proper- 
ties of the band structure leads to the formation of pro- 
tected metallic states at the surface. Soon after their 
theoretical prediction^!- as a generalization of Haldane's 
Chern insulator— their two-dimensional (2d) variant — 
the quantum spin Hall insulator — was realized in quan- 
tum well heterostructures by Konig et alA. Later, three- 
dimensional (3d) varieties of this novel form of insulator 
were predicted^"— and again realized afterwards^. 3d 
TIs come as strong and weak TIs. Strong TIs (STIs) 
have a metallic surface which is protected against local- 
ization due to its helical nature. Band topology enforces 
that the surface metal can be described as an effective 
2d Dirac theory of massless electrons. Soon after the re- 
alization of 3d STIs it was understood that the surfaces 
of these materials could give way to an unconventional 
Hall response whenever the surface is gapped by pertur- 
bations which break time-reversal symmetry (TRS). This 
is related to the axion quantum electrodynamicsiSrJ^ and 
the variation of the 0-angle at an interface between a 
bulk material {9 ^ 0) and vacuum (0 = 0). STIs are 
characterized by the quantized value 9 = n. When a gap 
is opened, a dissipationless surface conductivity with a 
half- integer Hall conductivit y^^i^^ , a^y = arises. Im- 
portantly, this value is not quantized if both TRS and 
inversion symmetry (IS) are broken in the bulk. 

Recently, Briine et al? have experimentally observed 
a quantized Hall effect for films of HgTe which can be 
considered 3d STIs, with the magnetic field applied per- 
pendicular to the film. Both in the experiment and in two 
theoretical studies^"— a field configuration was used, 
where only two faces of the crystal were gapped. 

In this paper we investigate the Hall response of a 3d 
STI in the presence of a bulk magnetic field oriented such 
that all surface excitations are gapped. We first study 
how the 0-term characterizing the 3d bulk, the surface 
Hall conductivities, the edge channels, and the quantum 
Hall effect observed by contacting those edge channels are 
related to each other. We then study a concrete lattice 
model (inspired by the HgTe band structure) to investi- 
gate how the location and number of edge channels can 



be controlled. 

We consider a single crystal STI with planar sur- 
faces in the presence of a magnetic field in a configu- 
ration, where both the bulk and all surfaces — up to 
one-dimensional edge channels discussed below — are 
gapped. The macroscopic electrodynamic response is de- 
scribed by an effective topological field theory in (3-f 1) 
dimensions defined by the actioniiii^ii^iii 
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This term is referred to as the 0-term, well-known from 
the field theory of axion electrodynamics"'^". Several 



properties 
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are important in the following: (i) The 



value of 9 is defined modulo 27r. (ii) Se is an integral over 
a total derivative, i.e., it has no effect for 9 = const., but 
matters at interfaces and surfaces, where 9 changes (see 
below). (Hi) Within the bulk, 9 = 91, is constant (we de- 
fine < < 27r to avoid ambiguities). In the presence 
of either IS or TRS, 0b assumes the values or tt. As 
TRS band insulators have 9^, = and STIs have 9^ = tt, 
this can be used to classify STIs (even in the presence of 
interactions). In the absence of both IS and TRS, 9^ is 
not quantized, and a Tl can adiabatically be transformed 
into a band insulator. 

Varying Eq. ([1]) with respect to one obtains dissi- 
pationless currents 
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located at the surface across which 9 changes. For a ho- 
mogeneous electric field at a surface Si one obtains the 
surface Hall conductivity by integrating the current per- 
pendicular to Si. Here it is important to realize that in 
the vacuum 9 = 2?™ with integer n. We find it con- 
venient to associate with each surface Si an integer rii 
defined by the value of 9 obtained by approaching the 
surface (see below for a numerical calculation). Defining 
the Hall conductivity of Si by measuring the current 
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FIG. 1: (Color online) a,b) Two possible edge channel con- 
figurations which are projected onto a plane in c) and d), 
respectively. On each face the integer Ui is written which 
characterizes the surface and from which the number of edge 
channels can be obtained using Eq. ([5} (for the planar graphs, 
n; = on the outside). The two configurations can be ob- 
tained by small rotations of the field, see Fig.|4l e) Orientation 
of line integrals used in Eq. ([5]). f) Numbering of edges and 
coordinate system. 



in the direction fiiX E for an electric field parallel to the 
surface, one obtains from Eq. ([5]): 
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For = one recovers the integer quantization of a 
conventional quantum Hall system and for Oh — the 
well known shift obtained for the 2d Dirac equation. In 
general, any value of 9\y is allowed. For a 3d sample 
with N surfaces the N integers rii and the constant 9^ 
completely characterize the macroscopic Hall response. 

As a simple application we can consider a quasi-2d 
slab with surfaces 5*1 and 52, and an electric field which 
is the same on the top and bottom of the sample. As 
fii = —712, the total 2d Hall response is given by 
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yielding the well-known integer-quantized form which is 
independent of the non-quantized 0^,- 

This analysis does, however, not answer the question 
what is actually measured when metallic contacts are at- 
tached to a 3d sample to obtain the Hall response. Inter- 
estingly, Chu, Shi and SheiJ^ argue that in the presence 
of dissipative surfaces (not considered by us) it might 
be possible to observe directly an approximately half- 
quantized in a four-terminal measurement. But in 
general, neither nor ajf are measured. Instead, one 
has to investigate a network of gapless current-carrying 
Id channels (also bent, non-planar quantum Hall junc- 
tions show similar networksi^). Each chiral Id channel 
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carries exactly one conductance quantum, / = j^A^, 
where $ is the scalar potential. To find those, we use 
Eq. ^ to calculate the current through a 2d mani- 
fold whose boundary is given by an equipotential line, 
$(r) = $c- From Stokes' theorem one obtains assuming 



that 0{r) is singular at most at one point on the manifold 
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with — $s — $0, where <I>s is the potential at the sin- 
gularity. Therefore, the number of Id channels is given 
by the (integer) winding number W of the 0-term, and 
their location tracks the singularities in 6{r) associated 
with yy ^ 0. The line integral is taken along the bound- 
ary of the manifold in clockwise direction looking along 
the direction of / [see Fig.[l]e)]. For crystals with flat sur- 
faces, the channels are localized along the edges, and W 
is obtained from the integers characterizing the adjacent 
surfaces, W = rii — nj, where i (j) is the surface with 
surface normal parallel (antiparallel) to the integration 
direction, respectively. 

Note that the winding numbers are quantized for any 
value of ^b- This reflects that the quantization of the 
Hall effect (based on charge quantization) is a more ro- 
bust concept than the Z2 quantization used to classify 
TIs (which relies on TRS). Deforming the TI into a ring 
(with a macroscopic cross section) one can directly apply 
Laughlin's gauge arguments^ for an alternative proof of 
the quantization of the Hall effect measured by contact- 
ing those edge channels. 

To determine how different networks of edge channels 
can be realized, we study a concrete example. As a min- 
imal modcl^ for a TI we consider four bands in a cubic 
lattice using orbitals inspired by strained 3d HgTe- (dif- 
ferent orbitals are important for 2d HgTe) . We introduce 
four basis states^- within the unit cell: 

|i)-|i?i,t>-|.s,t>, |3> = |i?i,;) = |,s,;), 

|2) - \LHA) = i\p^+ipy.,i) - 2|z,t»/V6, 
|4) - \LH,i) = -Hp, - tpyA)+2\z,i))/V6. 

In this basis, TRS is implemented by Q = —iaz sqK, 
where az acts in spin space, and the 2x2 identity matrix 
So acts on the orbital degrees of freedom {K is complex 
conjugation, a and s are standard Pauli matrices). Parity 
is implemented via P — ao 'Si Sz- Two- fold rotation sym- 
metries along the x-, y-, and z-directions are generated 
by Rx{t^) = ifJx ® Sz, Ry{T^) = —i(^y ® Sz, and Rz{t^) = 
ioz ® So- The 4x4 Hamiltonian can conveniently be ex- 
pressed in terms of the identity, F*^, five Dirac matrices 
F'^, and their ten commutators F"'' = [F°, F'']/(2i) which 
satisfy the Clifford algebra, {F",F''} = 2(5a,6F°. Using 

p(l,2,3,4,5) = (1 ^ ^ ^ ^ g^) 

we find as a minimal model 
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where m is the tuning parameter and i = 1 is the hopping 
amplitude. Ai and A2 parametrize the breaking of IS, 
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FIG. 2: (Color online) Number of edge channels along edges 
1-6 defined in Fig. [T] e) for = 0.02, Ai = 0, B = 
0.001(1, cos sin (p)/\/2, as a function of the rotation angle 
ip obtained from Eqs. (f9l 1101 [TT|) . The sign indicates whether 
the current runs parallel or antiparallel to the x, y and z di- 
rections. For the two configuration indicated by gray vertical 
lines, the edge configuration is shown in Fig. [T] a) and b). 
The number of channels on the six other edges is obtained by 
inversion symmetry (assumed here). 
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FIG. 3: (Color online) a^y^i) and 0i = 2Tvj + ^ J2j<i '^^vU) 
as a function of distance from the surface for an infinite slab 
in the xy plane (20 layers, <l>i = 0, »+Bz = 0.25). The 
integer j is chosen such that < S < 27r in the center. 
While each surface conductivity and 6h are not quan- 
tized, ajf = T.i^^^vi'i') and, equivalently, W = 6^(20) - 61(0) 
are quantized. W jumps from 2-k to at Ai ~ g^Bz ~ 0.25, 
where mz in Eq. ^ changes sign and one surface undergoes 
a QH transition. 



Ai also breaks TRS. For systems without inversion sym- 
metry (as HgTe), A2 7^ and a finite Ai will be induced 
by magnetic fields. Orbital effects of the magnetic field 
are described by the Peierls factor exp{iAn.n+kj) defined 
on the link from site n to the neighboring site n + Sj. 
The Zccman effect is obtained from 
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^ Blig+T'' - .g-r24) + ^Kg+F^^ ^ ^-p^)] (7) 

where Bz = I^bB and g^ = {gsi ^9lh)/'^ are linear 
combinations of the g factors of the El and LH sub- 
bands. For our plots we use m = —2, A2 = (A2 7^ 
does not change our results qualitatively), g^ — g^ = 1, 
and for a magnetic field in the j/z-plane we use the Lan- 
dau gauge, i.e., An^n+e, = 27r($j,z„ - $2y„)5ej ,e, , whcrc 
$i is the flux per unit cell in direction i in units of the 
flux quantum. 

For Ai = A2 = and vanishing magnetic fleld, the 
model describes a trivial band insulator for \m\ > 3 
and a STI for < |m| < 3. At m = ±3 the bulk 
band structure is characterized by a 3d Dirac point. For 
< Am = m -I- 3 < 1 and A^, B^, < Am one can 
obtain analytically the theory for a surface with normal 
vector h. First, one obtains the surface state for van- 
ishing Ai, i?^, $i from the standard bound state solution 
of the 3d Dirac equationi^ with space dependent mass 
niD = —Am sign(n-r), where a negative m^i for n r > 
mimics the vacuum. From this one obtains the 2d Dirac 
equation describing the surface of a TI^^. To linear order, 
the effects of A^ and i?^ are just given by the projection 
of H on the surface state, while the orbital effects are 
covered by minimal substitution. We obtain 

i/^-^f ^VFa-{p-^^A)~aofi + mza^ (8) 
mz =h-{^g-B%,~g-Bl,g+BI,)-Ai (9) 



where er ~ ((T^,cr^), p — {px,Py) denote the momenta 
perpendicular to h (shifted by a Bz dependent con- 
stant), is the Pauli matrix parallel to n, and /i de- 
notes the chemical potential. A is the vector potential 
of the normal field component appropriately transferred 
into the new basis. Note the unexpected dependence of 
mz on the direction of the magnetic field which allows 
to control orbital and Zeeman effects independently. For 
example, mz can be finite even for a magnetic field per- 
pendicular to the surface. 

It is an elementary exercise'^ to compute the spectrum 
of the massive Dirac theory in a field 



Ej = sign(j)Wa;2|j| + Ej=o = sign(B_L)™z- (10) 



with integer j and uj^ — vpy^2\eB±\. B± — B h de- 
termines the orbital contribution of the magnetic field. 
As long as B is not almost parallel to the surface, mz is 
only important for j = as m\ ^ (hbB)^ <^ Vp\eB±\. 

As the levels (fTO|) do not cross, the results for the in- 
tegers rii characterizing the surfaces, do not depend di- 
rectly on mz and are obtained from the QHE for Dirac 
fermions'^ and Eq. ([3|): 

n=^' + ^' l^^l fori?,<^<i?,+i. (11) 

Using Eqs. (|9I101 [TT]) one can therefore calculate the rel- 
evant winding numbers which determine the position of 
the edge channels, see Eq. ([5]) and Fig. [1] Note that /i 
can be controlled independently on each surface by ap- 
propriate gates and might also depend on details of the 
surface chemistry. As an example, we show in Fig. [2]how 
the number of edge channels changes when a magnetic 
field is rotated around the x axis. This allows to realize 
a wide variety of channel configurations. 
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FIG. 4: (Color online) Band structure of a beam of dimen- 
sions (oo X 20 X 20) as function of the momentum along the 
beam ($y = $^ = 0.018, g' Bl = = 0.28, Ai = 0.1). 

Flat parts in the central region represent Landau levels while 
dispersive parts are one-dimensional edge channels. Lower 
panel: Position of the edge channels for four chemical poten- 
tials marked above. Red (blue) denotes edge channels running 
parallel (antiparallel) to the x direction. For /x « 0, when mz 
dominates, B points towards/away from the edge channels, 
while a more conventional perpendicular arrangement is ob- 
tained in the orbitally dominated regime, /is. The last figure 
shows three edge modes moving in x direction on one edge, 
while three —x modes are distributed over the other three 
edges. Also two counterpropagating channels at the same 
edge can be obtained, e.g., for jl (not shown). 



The low-energy theory ([8|) with infinite cutoff does not 
account for the change of the 6'-term in the bulk, while 
it can be computed numerically by considering the full 
model ([B])- For an infinite slab in the a:y-plane, we have 
calculated from the Kubo formula the matrix of conduc- 
tivities, (Jxyih j), describing the current in the x direction 
in layer i when an electric field is applied in layer j in 
the y direction. For simplicity, we considered here only 
a Zeeman magnetic field. axy{i,j) is to a good approx- 
imation diagonal (as implicitly assumed in Eq. [1]), and 
we plot in Fig. [3] (Ta;j,(«) = J2j '^xyih j), i-c, the current 
distribution for a uniform electric field, and also the cor- 
responding layer-resolved 9i obtained from Eq. For 
finite Ai such that TRS and IS are broken, and the 
surface conductivities are not quantized, but the total 2d 
Hall response remains quantized as in Eq. (U). 

To illustrate the localization of the edge channels along 
the edges of the 3d crystal and to show how our find- 
ings relate to the band structure, we have calculated the 



properties of an infinite beam of dimension oo x 20 x 20 
using exact diagonalization. In the appropriate Landau 
gauge, the momentum along the beam is conserved and 
one obtains the band structure of Fig. 2] a) . To obtain 
the position of the conducting channels we plot in Fig. |4] 
the corresponding probability distribution: as expected 
they are located at the edges of the sample. By changing 
the chemical potential a sequence of quantum Hall tran- 
sitions is obtained. The first three figures show how the 
edge channels jump from one edge to the other, conserv- 
ing the number of channels. For other cases, the number 
of channels can also be increased as in a conventional 
quantum Hall transition and one can obtain more com- 
plex situations, see last panel. 

For a given edge channel configuration one can deter- 
mine the conductance measured when contacting those 
channels using (i) Eq. (O, (ii) Kirchhoff's law (current 
conservation at each corner) and (iii) that the electro- 
chemical potential remains constant when channels split. 
For the circuit shown in Fig.[l]d) one obtains $a = ^& = 

$d, = = - $a - ^I, $c = $a - 2^/, 

and $/ — — Defining the cross-conductances 

Gap = I/i'^p — this implies that not only, e.g., 

Ggb = /h and Gfa = 2e^//i can be measured but also 
Gcf = |e^//i. At the points where two channels meet, 
the power Pj is dissipated, with a quantized ohmic con- 
ductance l"^ /Pi which takes the values 2e^//i and ^e'^/h 
at the 'merging' corners above contacts e and f, respec- 
tively. Such a quantized heat source might have interest- 
ing applications for the study of heat transport in TIs. 

Can a topological insulator uniquely be identified by 
measuring the quantum Hall effect? For a fixed magnetic 
field this is not possible since there are situations where, 
for example, all surfaces can be gapped and there is not 
a single edge channel anywhere. The changing pattern 
of edge channels when field strength or direction or gate 
voltages are varied, can, however, serve as a fingerprint 
of topological insulators and their Dirac metal. The net- 
work of edge channels also opens new experimental op- 
portunities to study quantum Hall physics by separately 
addressing edge channels. Another interesting problem is 
the physics of fractional quantum Hall states, especially 
in combination with a non-trivial topology of surfaces 
and edges obtained by drilling holes into a 3d topological 
insulator. 
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